Hiding relativistic degrees of freedom in the early universe 
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We quantify the extent to which extra relativistic energy density can be concealed by a neutrino 
asymmetry without conflicting with the baryon asymmetry measured by the Wilkinson Microwave 
Anisotropy Probe (WMAP). In the presence of a large electron neutrino asymmetry, slightly more 
than seven effective neutrinos are allowed by Big Bang Nucleosynthesis (BBN) and WMAP at 2a. 
The same electron neutrino degeneracy that reconciles the BBN prediction for the primordial helium 
abundance with the observationally inferred value also reconciles the LSND neutrino with BBN by 
suppressing its thermalization prior to BBN. 



I. INTRODUCTION 

If the three known neutrinos (y e , v^, v T ) mix only with 
each other, all active neutrino degeneracies will equili- 
brate to close to the electron neutrino degeneracy via 
neutrino oscillations before BBN begins [1]. There are 
no caveats to this statement within the standard cos- 
mological model for the oscillation parameters relevant 
to atmospheric neutrino data and the Large Mixing An- 
gle (LMA) solution to the solar neutrino problem (which 
has recently been confirmed by the KamLAND experi- 
ment [2]). Thus, the magnitude of the electron neutrino 
degeneracy allowed by BBN is of special interest to any 
determination of, or constraints on, a lepton asymmetry 
in the universe. 

The authors of Ref. [1] emphasized that the degenerate 
BBN scenario, in which a small v e asymmetry 1 £ e ~ 0.2 
(which strongly affects the neutron to proton ratio and, 
as a result, the primordial 4 He abundance) is compen- 
sated by a much larger or v T asymmetry |£ p , r ^2 — 3 
(which affects the expansion rate) is excluded by this 
equilibration, leading to a limit £j| < 0.1 on all of the 
asymmetries. 

However, larger asymmetries may still be allowed if 
the effect of £ e on BBN is compensated by a new source 
of energy density other than that from the asymmetries 
themselves 2 . In particular, if there are additional degrees 



1 For a neutrino flavor a, an asymmetry L a between the numbers 
of v a and v a ("neutrino degeneracy") can be described by the 
neutrino chemical potential fi a or by the dimensionless degener- 
acy parameter £ a = fi a /T: 
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of freedom which do not mix with the three active neu- 
trinos, but which do contribute to the relativistic energy 
density at the BBN epoch, their presence can be hidden 
by an electron neutrino degeneracy. In this paper we 
quantify and explore the extent to which this statement 
is true. 

As an example, light Dirac neutrinos involve three 
additional S[/(2)-singlet neutrinos vr, the right-handed 
partners of the active states. Within the standard model 
(extended to include the right-handed neutrinos) these 
states have no interactions other than those associated 
with neutrino mass, i.e., they can only be produced by 
mass effects of amplitude m v /E or the associated tiny 
Yukawa couplings. These are too small to yield signifi- 
cant number densities at the BBN epoch for masses in 
the sub-eV range. However, light Dirac neutrinos are 
especially motivated in models involving new gauge in- 
teractions beyond the standard model (which may forbid 
the traditional seesaw mechanism) [3] . In that case they 
need not be sterile with respect to the new interactions, 
and may be efficiently produced prior to BBN [4]. This 
was recently studied in detail for a class of models with 
Z' couplings motivated by Eq [5]. The production and 
decoupling of the right-handed neutrinos, and the sub- 
sequent dilution of their number density by annihilation 
of decoupled heavy particles and color confinement fol- 
lowing the quark-hadron phase transition was calculated. 
The results are dependent on the specific Z' couplings 
and mass but, typically, the energy density at BBN as- 
sociated with the vr could be equivalent to that carried 
by 3, 2, or 1 active neutrinos for a Z' mass of 500 GeV, 
1 TeV, or 2 TeV, respectively, with the constraints be- 
coming weak or disappearing for the specific couplings 
for which the Z' decouples from the vr. 

If sterile neutrinos that mix with the active ones exist, 



Note that L a « 0.684£ Q for £ Q < 1. 
2 The faster expansion rate due to the new energy density could 
in principle also affect when or whether the equilibration occurs. 



However, the considerations in this paper are independent of 
equilibration. 
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whether or not flavor equilibration occurs is not known. 
We shall only consider this case in the context of the 
LSND [6] neutrino for which a large electron neutrino 
degeneracy suppresses its mixing with the active neutri- 
nos [7]. There has been considerable interest in sterile 
neutrinos as a way to account for the antimuon neutrino 
to antielectron neutrino oscillation signal observed in 
the LSND experiment corresponding to a neutrino mass- 
squared difference scale of about O.f — 1 eV 2 [8]. With 
three neutrinos there are only two independent mass- 
squared difference scales, that are fixed by the observed 
atmospheric and solar oscillations to be about 2.5 x 1 0~ 3 
eV 2 [9] and 7xfCT 5 eV 2 [10], respectively. Consequently, 
oscillations to a sterile neutrino have been invoked to 
explain the higher mass-squared difference scale of the 
LSND effect. The MiniBooNE experiment [11] now run- 
ning at Fermilab is designed to confirm or refute the 
LSND effect. Analyses with four neutrinos do not provide 
a good global fit to the data [12] , but the possibility exists 
that the inclusion of certain small parameters, hitherto 
neglected in the analysis, may improve the fit [13]. 

We account for additional relativistic (at BBN) degrees 
of freedom ll X v by normalizing their contribution to the 
energy density to that of an equivalent, non-degenerate, 
neutrino flavor [14], 

Px = AN vPv = -AN v p^. (2) 

o 

The effective number of neutrinos is N v = 3 + AN^ . 

Throughout this work the following primordial values 
are adopted for the 4 He mass fraction Y and for the deu- 
terium abundance j/d = 10 5 (D/H): 

Y = 0.238 ± 0.005, t/ D = 2.6±0.4. (3) 

For a detailed discussion of these choices, see Ref. [15]. 

II. THE ROLE OF £ e IN BBN 

The baryon-to-photon ratio 77 provides a dimension- 
less measure of the universal baryon asymmetry which is 
very small (rj = ub/u^ < 10~ 9 ) 3 . By charge neutrality 
the asymmetry in the charged leptons must also be of this 
order. However, there are no observational constraints, 
save those to be discussed here (see [16, 17] and further 
references therein), on the magnitude of any asymmetry 
among the neutral leptons (neutrinos) . A relatively small 
asymmetry between electron type neutrinos and antineu- 
trinos (|£ e | > 10 ~ 2 ), but large compared to the baryon 
asymmetry, can have a significant impact on the early 
universe ratio of neutrons to protons, thereby affecting 
the yields of the light nuclides formed during BBN. The 



3 We define r?io = 10 10 r/, to facilitate the use of numbers of order 
unity. 
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FIG. 1: Isoabundance curves for D and 4 He in the £ e - f?io 
plane for N v = 3. The nearly horizontal curves are for 4 He 
(from top to bottom: Y = 0.23, 0.24, 0.25). The nearly ver- 
tical curves are for D (from left to right: j/d = 10 5 (D/H) = 
3.0, 2.5, 2.0). The data point with error bars corresponds to 
j/d = 2.6 ± 0.4 and Y = 0.238 ± 0.005. 



strongest effect is on the 4 He abundance, which, during 
BBN, is neutron limited. For £ e > 0, there is an ex- 
cess of neutrinos (v e ) over antineutrinos (P e )i and the 
two body reactions regulating the neutron to proton ra- 
tio (n + v e <-> p + e~ , p + D e «-> n + e + ), drive down the 
neutron abundance; vice- versa for £ e < 0. The effect of 
a non-zero v e asymmetry on the relic abundances of the 
other light nuclides is much weaker. This is illustrated in 
Fig. 1, which shows the D and 4 He isoabundance curves 
in the £ e — 7710 plane for N v = 3. The nearly horizontal 
curves reflect the weak dependence of Y on the baryon 
density, along with its significant dependence on the neu- 
trino asymmetry. In contrast, the nearly vertical curves 
reveal the strong dependence of j/d on the baryon den- 
sity and its weak dependence on any neutrino asymmetry. 
This complementarity between j/d and Y permits the pair 
{77, AN V } or the pair {77, £e} to be determined once the 
primordial abundances of D and 4 He are inferred from 
the appropriate observational data. 



III. ANALYSIS 

In the "standard" case of N v = 3, the observed abun- 
dances of D and 4 He may be used to identify the al- 
lowed region in the rj — £ e plane. In addition, the cos- 
mic background radiation (CBR) data from WMAP [18], 
provide complimentary information for the 77 distribu- 
tion (for different N v [15]), which will further restrict 
the allowed 77 range. In Fig. 2 are shown the lcr and 
2a contours in the 7710 — £ e plane for N v = 3 from a 
joint CBR (WMAP) - BBN fit using the adopted D and 
4 He abundances. Although very large when compared to 
the baryon asymmetry, the allowed values of £ e are suffl- 
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ciently small (|£ e | < 0.1) that the "extra" energy density 
contained in such neutrinos is negligible 4 , justifying the 
claim that these results do, in fact, represent N v = 3. 
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FIG. 2: The la and 2a contours in the rjio — £ e plane for 
N v = 3 from a joint CBR - BBN fit using WMAP data and 
the adopted D and 4 He abundances. The cross marks the 
best-fit point (r?io,£ e ) = (6.16,0.044). 



A. The LSND neutrino 

If the LSND results [6] are to be accounted for by 
the mixing of a "sterile" neutrino with one of the active 
neutrinos, it is of importance to consider whether the 
LSND neutrino will be populated, through mixing with 
the active neutrinos, sufficiently early in the evolution of 
the universe that it contributes AY„ = 1 to the rela- 
tivistic energy density prior to BBN. Given the parame- 
ters required by LSND, in the absence of any significant 
asymmetry among the active neutrinos the LSND neu- 
trino would, indeed, have come into thermal equilibrium 
sufficiently early [7]. In this case the new "standard" 
model would have N v = 4. However, in the absence of 
a significant neutrino degeneracy, N v — 4 is strongly ex- 
cluded [15]. 

In contrast, as the asymmetry in the active neutri- 
nos increases, the sterile-active mixing is delayed [7]. For 
£ e ^ 0.01 — 0.1, mixing occurs after the active neutrinos 
have decoupled from the e^ - 7 plasma. In this case, al- 
though the sterile neutrino does still mix with the active 
ones, the energy is shared among all the neutrinos and 
N v remains equal to 3. It is interesting, even amusing, 
that the range of £ e required by BBN for N v — 3 is, in 
fact, in excellent agreement with that required to delay 
the thermalization of the LSND neutrino (£ e > 0.01 - 



4 ajv„(0 = J2i [f (#) 2 + f(#) 4 ] ~ °- 013 ; scc > e -f- [is]. 



0.1), thereby keeping N v ~ 3. Thus, the same values 
of £ e that reconcile the low value of Y with the BBN 
prediction (for N v = 3) also reconcile the LSND sterile 
neutrino with BBN by delaying its mixing with the active 
neutrinos. 



B. Arbitrary N v 





1 1 1 1 




1 | 
o.U 






/7.0 






/ 3.0 






' / 5.0 

1 







1 



-2 2 4 



FIG. 3: The approximate £ e - AY„ relation corresponding 
to yr> — 2.6 and Y = 0.238. The numbers shown are the 
corresponding values of 7710 at those points. 



Let us set aside the LSND neutrino, assume that Y„ = 
4, and ask if BBN and the adopted primordial abun- 
dances can be reconciled (along with the CBR) through 
the appropriate choice of a nonzero asymmetry in the 
active neutrinos. 

From BBN with fixed abundances of D and 4 He, as 
N v increases, so too will the best fit values of rj and £ e . 
It is straightforward to understand this effect. As N v in- 
creases, the early universe expands more rapidly, leaving 
less time to burn D. As a result, for a fixed baryon to pho- 
ton ratio, j/d increases with increasing N u 5 . To reduce 
?/d back towards its observed value requires increasing the 
baryon density. The combination of an increased baryon 
density and N v — 4, raises the predicted primordial abun- 
dance of 4 He, requiring a larger £ e to reconcile the BBN 



J This effect on the primordial deuterium abundance is often, but 
not always [15, 19], ignored. However, for |AJVi/| ^ 1, it is not 
negligible. 
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predictions with the data. Indeed, it is clear from Fig. 1 
of Kneller et al. [17] that for any (reasonable) choice of 
N v there is a pair of {77*0, £*} values for which perfect 
agreement with the observed D and 4 He abundances can 
be obtained (xbbn = 0)- An approximation to this £ e 
- AN„ relation [20] is shown in Fig. 3 for j/d = 2.6 and 
Y = 0.238; the numerical values along the curve are the 
corresponding values of 7710- However, this correlation of 
increasing 77 with increasing N u is broken by the inclusion 
of WMAP data through the CBR-imposed constraint on 
the rj - N v relation [15]. As a result, while for N v = 4 the 
allowed values of £ e increase from that for N v — 3, the 
shift in rj is much more limited. In Fig. 4 is shown the 
analogue of Fig. 2 (the confidence contours in the rjw — £ e 
plane) for N v = 4. 
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FIG. 4: Same as Fig. 2, except N u — 4. The cross marks the 
best-fit point (7710,6s) = (6.22,0.096). 



It is clear that from BBN alone, reconciling the primor- 
dial D and 4 He abundances with the BBN predictions for 
large N„ would require too large a baryon to photon ra- 
tio to be consistent with the rj - N v constraints from 
WMAP. The parameter pair {AN v ,(, e } can therefore be 
constrained using the {AiV„,77} likelihood distribution 
obtained [15] from the WMAP data [18] along with the 
{£ e , AN„,7j} likelihood distribution from our BBN analy- 
sis, by forming the joint likelihood and marginalizing over 
rj with the prior 71(77). Adopting a flat p.d.f. for 77(77) the 
{AiV„,£ e } constraints are derived after calculating (up 
to a multiplicative constant), 

L(£ e ,AN„) = J L BBN (t e ,AN„,r))xLwM AP (AN v ,r))dri. 

(4) 

This analysis results in the allowed regions in the AN U - 
£ e plane shown in Fig. 5, providing limits on the values 
of AiV„ that can be accommodated by non-zero £ e with- 
out conflicting with the 77 and AN V values allowed by 
WMAP. 
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FIG. 5: The la and 2a contours in the AiV„ - £ e plane 
using WMAP data and BBN with the adopted D and 4 He 
abundances. The cross marks the best-fit point (AiV„,£ e ) = 
(0.22,0.057). 



IV. CONCLUSIONS 

Big Bang Nucleosynthesis places strong constraints on 
any new source of relativistic energy density present dur- 
ing the BBN epoch. Examples include the energy density 
associated with large asymmetries |^, T | between ji or t 
neutrinos and their antincutrinos, with the right-handed 
components vr of light Dirac neutrinos if they can be 
efficiently produced by new interactions, or with sterile 
neutrinos which mix with active neutrinos, as suggested 
by the LSND experiment. In the absence of a significant 
neutrino asymmetry the data constrain the equivalent 
number of new neutrino flavors AN V to AiV^ < 0.1 — 0.3 
(assuming that AiV„ > 0) [15]. 

BBN also constrains the electron neutrino asymmetry, 
£ e , which directly affects the neutron to proton ratio and 
therefore the 4 He abundance. For AN V = we find £ e < 
0.1, as shown in Fig. 2. However, there remains the possi- 
bility of compensation between the effects of AjV,, and £ e , 
which can relax the constraints on both. The expanded 
ranges, consistent with WMAP and BBN, are shown in 
Fig. 5. For £ e in the range —0.1 < £ e < 0.3, values of 
AN V in the range —2 < AN U < 5 are permitted. After 
marginalizing over £ e , the 2a range for AN„ (obtained by 
identifying those values of AN U for which the likelihood 
is above L max /e 2 ) is -1.7 < AN U < 4.1. 

The possibility of compensation between a small £ e 
and a large |^ )T | is now excluded given the oscillation 
parameters inferred for the atmospheric neutrinos and for 
the LMA solution for the solar neutrinos because of the 
equilibration of the three £i [1]. However, compensation 
between £ e and other types of relativistic energy density, 
such as neutrinos which do not mix with active ones, 
remains a possibility. We have explored this possibility 
in detail in this paper finding, for example, that AN V = 1 
is consistent with the BBN and the adopted primordial 
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abundances provided that £ e ~ 0.1, as shown in Fig. 4. 
Another consequence of our analysis is that even three, 
light, Dirac neutrinos (AY„ = 3) can be accomodated 
provided that £ e ~ 0.2; see Fig. 5. 

Mixing between sterile and active neutrinos must 
be considered separately, because £ e can shift the 
time/temperature at which the sterile neutrinos are pro- 
duced [7]. Values of £ e in the range 0.01 — 0.1 (which are 
allowed for AiV„ <~ 0; see Fig. 2), are sufficient to keep 
N v = 3 and thus reconcile the LSND neutrino with BBN. 

The compensation between A7V„ and £ e explored here 
is, perhaps, fine-tuned. It should not be dismissed, how- 
ever, because of the significance of these constraints for 
particle physics (in particular, for models of neutrino 
mass and mixing), and also because of the importance 
of AN V for the subsequent evolution of structure in the 
universe, as probed by the CMB and large scale struc- 
ture surveys. Also, a large £ e of the magnitude relevant 
for BBN (e.g., £ e ~ few x 10 _1 ) is huge compared to 
the baryon asymmetry, so even if it is unlikely, a value in 
this range would be of profound significance for particle 



physics and cosmology. While in many models of lepto- 
genesis the asymmetry among the neutral leptons is tied 
to that in the charged leptons and, hence, is very small 
(£e ~ V ~ 10~ 9 ), there are mechanisms for generating 
the large asymmetry required here [21]. 
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